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The phase transition between charge- and spin-density-wave (CDW, SDW) phases is studied in the 
one-dimensional extended Hubbard model at half-filling. We discuss whether the transition can 
be described by the Gaussian and the spin-gap transitions under charge-spin separation, or by a 
direct CDW-SDW transition. We determine these phase boundaries by level crossings of excitation 
spectra which are identified according to discrete symmetries of wave functions. We conclude that 
the Gaussian and the spin-gap transitions take place separately from weak- to intermediate-coupling 
region. This means that the third phase exists between the CDW and the SDW states. Our 
results are also consistent with those of the strong-coupling perturbative expansion and of the direct 
evaluation of order parameters. 
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It has been pointed out that the phase transition 
between the charge-density-wave (CDW) and the spin- 
density-wave (SDW) phases of the one-dimensional (ID) 
extended Hubbard model (EHM) has curious properties. 
The Hamiltonian of this model is given by 

W = -*I>U+i, s +H.c.) 

is 

+U n^riji + V TtjTit+i ■ (1) 



Luttinger (TL) liquids |16|. In the TL liquid (bosoniza- 
tion) theory, the continuous fermion fields are defined by 

Cis -> 4>l,s(x) + ipn, s (x) with 



lp r s (x) = = c irk F x c i/\/2-\r(<l) +s(l>„)-8 D -s8„} ^ 



(2) 



From the strong coupling theory, it has been shown that 
the first-order transition takes place near U = 2V MSJ^. 
On the other hand, the weak coupling theory predicts 
that the transition is the second order on the same line 
This means that there exist a crossover between 
these two transitions while, in more than two dimension, 
the transition is the first order for any strength of the in- 
teractions. Many authors investigated this transition by 
analytical |(| and numerical approaches [||,(7|j|,||,[l(J . Nu- 
merical analysis has been done by the direct estimation of 
the CDW order parameter and by the real-space 

renormalization group technique 0. However, due to 
the lack of precision of the analysis, the property of this 
transition is still left to be ambiguous. 

In this Letter, we study this transition by investigat- 
ing level-crossings of excitation spectra 1 1 1 iajl3,n4,|5[ 



where r = R, L and s =|,| refer to +, — in that or- 
der, is the Fermi wave number. The field <p v and 
dual field 9 U of the charge (y — p) and spin (y = a) 
degrees of freedoms satisfy the relation [4> lJl {x),9 u (x')\ = 
— m 5^ sign (a; — x')/2. Then the effective Hamiltonian 
for the system with length L is given by the sine-Gordon 
model: 

U = E ?T I dx [ K v(W 2 + K-\d x ct>u) 2 ] (3) 
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which enable us to determine the phase boundary with 
high accuracy from the numerical data of finite-size clus- 
ters. Using this technique, we determine the phase 
boundaries as shown in Fig^. Then we conclude that 
two transitions of charge and spin degrees of freedom oc- 
cur independently, and the third phase exists between 
the CDW and the SDW phases. 

First, we briefly explain the level-crossing method used 
in our analysis. In general, the low-energy behavior 
of ID electron systems can be described as Tomonaga- 



where v v and K v denote the sound velocity and the Gaus- 
sian coupling, respectively, for each sector. The non- 
linear terms in eq.(|^) are originated from the Umklapp 
(charge) and the backward (spin) scattering effects which 
can cause the charge- and the spin-gap instabilities. If 
these non-linear terms are irrelevant, the excitation spec- 
tra and their wave numbers in this system are given by 
the relation 
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k = -j~( s p + s <?) + 2m p k F , 



(4) 
(5) 



where x v = [n1/K v + mlK l/ )/2 are the scaling dimen- 
sions and s u = n u m u are the conformal spins. The in- 
tegers n v and m„ are quantum numbers for the particle 
numbers and the current excitations, respectively. The 
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scaling dimensions are related to the critical exponents 
for the correlation functions as 



(Oi(r)Oi(r')) ~ | r -/|-2(^+^). 



(6) 



Therefore, there is one to one correspondence between 
the excitation spectra and the operators. 

In our analysis we turn our attention on excitation 
spectra which correspond to the following operators: 



O vl = \/2cos(V2^), 
O v3 = exp(±iv / 2^). 



(7a) 
(7b) 
(7c) 



For the spin sector (y = a) which have an SU(2) sym- 
metry, the operators (f7a| ) and j7b| , [7^ ) form singlet and 
triplet states, respectively. In this case, the level cross- 
ing of these spectra (x a i = x a 2,x a3 ) gives the spin-gap 
phase boundary On the other hand, the charge 

sector (v — p) is U(l) symmetric. Then the opera- 
tors of eqs.(0) correspond to "Neel", "dimer", and "dou- 
blet" states in that order, borrowing the terminology of 
anisotropic spin chains. In this case the level-crossing 
of "Neel" and "dimer" excitations {x p \ = x P 2) gives the 
Gaussian transition OJl7j], which means a second order 
transition between two massive states with different fixed 
points (g 3 ± — > ±oo). Note that these excitation spectra 
can be extracted when we choose anti-periodic boundary 
conditions (BC = — 1, see Tab. Jh, reflecting the selection 
rule of the quantum numbers E3 . In the weak-coupling 
limit, the Gaussian and the spin-gap transition take place 
simultaneously, because gi± = g 3 ± = U — 2V. However, 
there is no guarantee for the synchronization except for 
this limit. 

In addition to the above effective Hamiltonian (^), 
there exists the following Umklapp operator transferring 
finite spin 
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dx cos 



(8) 



In the weak coupling limit, the coupling constant is as- 
signed as <73|i = —2V, so that it remains finite on the 
U = 2V line. Therefore, we should consider the possibil- 
ity that the charge and the spin degrees of freedom is not 
separated, and a direct transition between the CDW and 
the SDW phases takes place. To examine this possibility, 
we also observe the level-crossing of excitation spectra of 
the CDW and the SDW operators which consist of both 
charge and spin components (see Tab. |). These ex- 
citation spectra can be obtained in periodic boundary 
conditions (BC = 1) with the wave number 2/cp = 71 "■ 

The excitation spectra correspond to the above opera- 
tors can be identified according to their discrete sym- 
metries. The wave functions for the excited states 
can change their signs under particle- hole (C: a s <-» 



(— iyc\ s ), space-inversion (V: R<->L), and spin-reversal 
(T: T^i) transformations. It follows from eq.(||), that 
the phase fields <\> v change by these transformations as 
follows: 



C,V 
T 



(9a) 
(9b) 



Here, CP = 1 is always satisfied, so that independent dis- 
crete symmetries are V and T. The boson representation 
of the operators and their symmetries are summarized in 
Tab. | ]l9(| . In the present numerical calculation based 
on the Lanczos algorism, the identification is done by 
projecting the initial vector as 



|* tai t) = 5(l±P)(l±T)|<) 



(10) 



where the signs in front of the operators correspond to 
their eigen values, and \i) is some configuration which 
satisfies V,T\i) ^ ±|*)- Furthermore, \i) is classified by 
the wave numbers k — 0, tt. 

The critical lines obtained by the above explained way 
with the exact diagonalization of the L = 8, 10, 12, 14 
systems are shown in Fig. ^. For the Gaussian transi- 
tion line, the finite-size effect is small for all region (see 
Fig.^(a)). On the other hand, for the spin-gap transition 
line, the finite-size effect is small in the weak-coupling 
region (see Fig.|^(b)), but large in the intermediate- and 
the strong-coupling regime. The CDW-SDW transition 
line lies between the above two lines [^o|. Its finite-size 
effect is large for all region. 

In order to check the consistency of our argument, we 
confirm the relations between the scaling dimensions for 
each instability. The relation on the Gaussian critical 
line |l^] and that near the spin-gap transition Jl^,|ljl are 
given by 



Xpl + X p2 

2— Xp3 

X al + 3x a2 ,3 



(11) 

(12) 
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The numerical results are shown in FigsJ^ and |]. The 
Gaussian and the spin-gap transition lines satisfy the 
consistency of our theoretical scheme from the weak- to 
the intermediate coupling region. 

Moreover, in order to back up the present result from 
the strong coupling theory, we compare above result to 
the critical line obtained by the strong coupling expan- 
sions. In this limit, the phase boundary between the 
CDW and the SDW states can be determined by equat- 
ing energies of these states. This calculation has already 
been done by Hirsch Q and van Dongen p] up to second 
and forth order, respectively, using the Bethe-ansatz re- 
sults. Among the three transition lines we assumed, van 
Dongen's result shows good agreement with the Gaussian 
transition in the charge part up to U/t ~ 6. We should 
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also note that our Gaussian critical point agrees with the 
Cannon et aZ.'s result obtained by the direct evaluation of 
the CDW order parameter: V c /t = 1.65±g*^ for U/t = 3 
V c /t = 2.92 ±0.04 for U/t = 5.5 [| (see Fig.' | and |(a)). 

From the above evidence, we conclude that the actual 
transition near U = 2V line is not direct transition be- 
tween the CDW and the SDW states, but independent 
Gaussian and spin-gap transitions at least from the weak- 
to the intermediate-coupling region. In the strong cou- 
pling regime, these two boundaries approach and coincide 
at the finite strength of the coupling. Unfortunately, in 
the present analysis, we can not determine this point, 
but it is considered to be identical to the crossover point 
between the second and the first order transitions. In 
this way, our analysis suggests that the crossover along 
the U — 2V line is closely related to the validity of the 
charge-spin separation. 

Our result also means that there is a finite region with 
charge- and spin-gapped state which has different sym- 
metries from the CDW state. This third phase is consid- 
ered as a bond-order-wave (BOW) state with LRO which 
is characterized by the following operator: 



Cbow 



(-1)' 



( C j+l,s 



(13) 



Therefore the direct evaluation of this order parame- 
ter may possible. The existence of the BOW state can 
be more clarified by extending the EHM. According to 
the bosonization analysis of the EHM with correlated- 
hopping interactions a finite BOW region remains 
even in the weak coupling limit (gi± < 53_l). On the 
other hand, if the Gaussian and the spin-gap transitions 
take place in the opposite order (g±± > 53_l), there ap- 
pears a bond-spin-density-wave (BSDW) phase which has 
massive charge sector and massless spin sector. 

Finally, we refer to other examples of the crossover 
similar to that in the EHM. This type of phenomenon has 
also been observed in the transition between the singlet 
and the Haldane phases in the 5 = 1/2 frustrated spin- 
ladder model J22) . Our result would also shed light on 
such a transition. 

In summary, we have studied the CDW-SDW transi- 
tion in the EHM along the U — 2V line, and shown that 
there exist the Gaussian and the spin-gap transitions in 
the charge and the spin degrees of freedom, respectively, 
and there is the BOW state between them. The crossover 
from the second to the first order transition is suggested 
to be related with the validity of the charge-spin separa- 
tion. 
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FIG. 1. Possible three transitions (Gaussian, spin-gap, and 
CDW-SDW transitions) along U = 2V line of the EHM cal- 
culated in L = 8, 10, 12, 14 systems. The result of the strong 
coupling expansion agrees with the Gaussian transition. This 
means that the actual transitions are the Gaussian and the 
spin-gap transitions, and a BOW state exists between them. 
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FIG. 3. Product of the scaling dimensions on the Gaussian 
critical line (eq.([ll])). The TL liquid theory predicts the value 
takes 1/4. 
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FIG. 2. Size dependence of the critical points for the (a) 
Gaussian transition at U/t = 3, and the (b) spin-gap transi- 
tion at U/t — 2. The former agrees with the Cannon et al.'s 
result that V c /t = 1.65^° ^ §. 
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TABLE I. Discrete symmetries of wave functions which 
correspond to several excitation spectra (C: charge conjuga- 
tion, V: space inversion, T: spin reversal, and k: wave num- 
ber). The upper (lower) sign of boundary conditions (BC) 
denotes N/2 —odd (even) cases. The upper 12 states are 
"physical" states which appear in the same boundary condi- 
tions as the ground state. The lower 6 states are the "ar- 
tificial" ones which are extracted by twisting the boundary 
conditions respect to the ground state. 
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